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ABSTRACT 

A high-resolution, multi-level, primitive equation ocean model is used to examine the response 
of an idealized, flat-bottomed, eastern boundary oceanic regime on a beta-plane to both steady 
and seasonally-varying climatological wind forcing. The focus of the study is the California 
Current System along the coastal region, from 35° N to 45° N, off the Western United States. 
With steady equatorward wind forcing, a surface equatorward current and poleward undercurrent 
develop. Eddies form around days 60 and 70 with initial development in the northern region of 
the domain. The strong meandering current continues to grow throughout the 360 days of model 
time and can produce eddies that have wavelengths up to ~ 200 km and can propagate at least 
~200 km offshore. When the alongshore component of the temporally averaged seasonally- 
varying climatological wind forcing is used, there is a weak poleward undercurrent and 
equatorward surface current. There ts weak upwelling and very little eddy activity with the eddies 
only propagating to ~100 km offshore. When alongshore component of the time-dependent wind 
forcing with spatial variability in latitude 1s used, a surface equatorward jet, poleward 
undercurrent and eddies are generated. The eddies form throughout the domain in this 
experiment due to a competition between the B-plane effect and the continuous and stronger 
equatorward wind forcing in the southern portion of the domain. The eddies in this experiment 
propagate at least ~150 km from shore. When the full climatological winds are used, there is 
again a surface equatorward jet and poleward undercurrent generated. The surface currents and 
eddies reach speeds of ~65-85 cm s*!. The wavelengths of the eddies are ~300 km and can 
propagate to ~350 km offshore. These results are closest to the observed structure of the 
California Current System. Therefore, I conclude that both temporal and spatial variability in 
wind forcing are critical elements in the formation and maintenance of currents and eddies in the 


California Current System. 
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I. INTRODUCTION 

The California Current System is a complex system of eddies and meandemng jets 
superimposed on a mean equatorward flow extending to ~1000 km from the coast. There 
have been several major experiments including Coastal Upwelling Experiments (CUE-I 
and CUE-II), Coastal Ocean Dynamics Experiment (CODE), Ocean Prediction through 
Observations, Modeling and Analysis (OPTOMA), and Coastal Transition Zone (CTZ), 
designed to better understand this eastern boundary current. The purpose of this study is 
to explain the generation of the eddies and filaments found in the above experiments. 
Batteen e/ a/. (1989) demonstrated how a steady equatorward wind 1s able to generate a 
reasonable current, undercurrent and eddy field with a full primitive equation (PE) model. 
The process-oriented study described here extends the work of Batteen e/ a/. (1989) by 
using seasonal rather than steady winds to force the same PE model. These winds are 
then averaged into different components to understand the role these components play in 
the generation of mesoscale features. 

This thesis is organized as follows. Chapter II descnbes the region being modeled 
as well as the type of winds used to force the model. A brief look at some numerical 
models used in the past is at the end of this chapter. In Chapter III there is a descnption 
of the model and the design of the experiments. The results of the model experiments are 


presented in Chapter IV. A summary ts provided in Chapter V. 


Il. BACKGROUND 


A. REGIONAL DESCRIPTION 

Oceanic flow off California is classified as an eastern boundary current. It was 
described as a broad, weak, equatorward flow extending to ~1000 km from the coast 
(Sverdrup ef al., 1942). Source water is carried eastward in the West Wind Dnift of the 
North Pacific Gyre (Lynn and Simpson, 1987). This water forms a relatively cold, low 
salinity current along the west coast of North Amenca. The pnncipal equatorward flow is 
above 1000 meter depth with the bulk of the average 11.9 Sv transported equatorward 
occurring within 1000 km of the coast (Wooster and Reid, 1963). Numerous eddies and 
irregularities in the flow as well as the existence of coastal countercurrents have been 
observed (Wooster and Reid, 1963). Equatorward winds in the summer can cause 
equatorward currents and strong upwelling while a coastal, poleward countercurrent with 
no upwelling can be present in the winter (Reid et a/., 1958). Non-seasonal vanations in 
upwelling appear to be related to wind changes (Reid, 1960). 

The California Current System (CCS) can be broken into four general currents: the 
Californa Current (CC), the California Countercurrent (CCC), the California 
Undercurrent (CUC), and the Southern California Countercurrent (SCC). Since only the 
first three are included in the area of this modeling study, the SCC will not be further 
discussed. 

The equatorward CC has a core that 1s within 100-200 km of the coast and reaches 
a depth of nearly 300 m (Chelton, 1984). The average current speeds are less than 25 cm 
s-! (Reid and Schwartzlose, 1962). Core speeds often in excess of 75 cm s"! have been 


observed (Bnnk ef al., 1991). . 


The poleward CCC is observed inshore of the CC and is characterized by its 
seasonal variability (Chelton, 1984; Hickey, 1979). The flow of the CCC 1s equatorward 
from February through September and poleward from November through February 
(Hickey, 1979). When flowing poleward on the surface, it is called the Davidson Current 
(Chelton, 1984). Observations north of Cape Mendocino show that the flow can be either 
poleward or equatorward dunng the same month (Chelton er al., 1988). 

The CUC flows poleward over the continental shelf (Hickey, 1979; Reid, 1962). 
Its core has a mean velocity greater than 15 cm s~! with instantaneous speeds of up to 40 
cm s7! (Hickey, 1979; Reid, 1962). The CUC is 20 - 70 km wide and vertically extends ~ 
300 meters with its center 200 - 250 meters below the surface (Hickey, 1979; Reid, 1962). 

The colder temperatures found near the coast and observations of the CCC below 
200 m depth were originally attnbuted to a classical two-dimensional upwelling system 
closely related to the winds (Smith, R. L., 1968). This simple one-celled circulation of 
upwelling was later disputed by observations during CUE-] and CUE-2, where the 
vertical distributions of temperature, chlorophyll and optical parameters contradicted the 
simple two-dimensional theory (Huyer, 1983). 

The CODE program in 1981-1982 was designed to study and identify those 
dynamical processes which govern the wind-driven motion of water over the continental 
shelf (Beardsley and Lentz, 1987). Kosro (1987) used the shipboard Doppler acoustic log 
(DAL) to show that near-surface currents often deviated substantially from classical two- 
dimensional wind-dnven upwelling. DAL results indicated a close association between the 
complex temperature structures in satellite imagery and the presence of vigorous current 
structures including squirts, eddies, and countercurrents (Kosro, 1987). The OPTOMA 
program further showed the persistence of the eddy field and its potential for drawing the 


cold waters offshore (Mooers and Robinson, 1984). The CTZ program was the next 


major experiment designed to help understand the complex CCS. A major portion of this 
program was dedicated to resolving what caused the filaments (cold surface features) 
found dunng CODE and OPTOMA (Bnnk and Cowles, 1991). Dnfters showed core 
speeds often in excess of 75 cm s“! in the generally southward meandering jet during July 
to December 1988 (Brink ef al., 1991). When the winds were poleward and the Davidson 
Current was present, there was little eddy activity and no coherent flow or filaments were 
observed (Kosro ef al., 1991). 

CTZ was designed to study filaments. Yet in the EBC region many large-scale 
eddies were present. The newest program, "Weakly Nonlinear Accelerated Research 
Initiative", sponsored by the Office of Naval Research (ONR), 1s designed to try to 
understand more about the eddies in the CCS. Of particular interest is the question of the 


role of wind forcing in generating eddies in the CCS. 


B. CLIMATOLOGICAL WINDS 

The atmospheric circulation off the west coast of North America 1s dominated by 
the North Pacific High, which varies seasonally in position and strength (Huyer, 1983). In 
February the High is centered near 28°N, 130°W, while in August it migrates to its most 
northward location of 38°N, 150°W (Huyer, 1983). Nelson (1977) showed that the 
strongest equatorward wind stress migrates from ~259N in the winter to ~38°N in the 
summer. This shifting of the High causes the winds to split at ~40°N near Cape 
Mendocino during the winter. The northward winds blow around the eastern side of the 
Aleutian Low, north of the cape, while southward winds blow around the North Pacific 
High, south of the cape. In the summer, the winds are predominately southward on both 


sides of the cape (Nelson, 1977). 


Smith (1968) showed that open ocean upwelling can occur if there is positive wind 
stress curl away from the boundary of the coast. Bakun and Nelson (1991) found that the 
large-scale wind stress curl in the CCS is predominately cyclonic near the coast and 
anticyclonic offshore. As a result, Ekman convergence dominates near the coast, while 
Ekman divergence dominates offshore. There is a transition zone located approximately 
200 - 300 km off the coast between the cyclonic and anticyclonic curl which 1s associated 
with the offshore maximum of alongshore wind stress (Bakun and Nelson, 1991). 

There are two local maximums in the wind stress curl. The first is yearlong and 
associated with the Southern California Bight. The second extends south from Cape 
Mendocino to past Monterey Bay and occurs during the fall and winter (Bakun and 


Nelson, 1991). 


C. NUMERICAL MODEL STUDIES 

Numencal models have been used to help unravel some of the mystenes 
concerning the causes and behavior of the phenomena observed in the CCS. They can 
give the solution to complicated, non-linear problems that are intractable with standard 
analytical methods (Allen, 1980). 

Two-dimensional models, which ignored alongshore vanability, were a_ first 
attempt at using the computer to generate results that were reasonable when compared to 
available observations. For example, the two-layered stratified models discussed by Allen 
(1980) with a flat bottom produced a coastal jet in the upper layer and a countercurrent in 
the lower layer. When a shelf and slope topography were added to the same model, 
results showed that the countercurrent strength would be significantly reduced, consistent 


with available observations. 


As computers have advanced in speed and power, three-dimensional models have 
replaced the two-dimensional models and can be run with results obtained in a reasonable 
amount of time. For example, for eastern boundary current regions, two-layer and 
continuously stratified three-dimensional models have been used by Suginohara (1977) to 
study the effects of alongshore vanations in wind stress. The results have been shown to 
be consistent with forced coastal trapped wave theories. The results also show a 
pycnocline upwelling near the coast and the formation of fronts, consistent with available 
upwelling observations (e.g., Curtin, 1979). 

None of the models mentioned above simulated eddies, jets or filaments. 
Haidvogei ef al. (1991) was able to produce these results but used no wind forcing. 
Batteen e/ a/. (1989) had the first eastern boundary current model to simulate eddies and 
jets with steady wind forcing, followed by McCreary e/ a/. (1991), who used the same 
type of wind forcing. 

In this study, the model used by Batteen ef a/. (1989) is extended to include 
seasonal rather than steady wind forcing. The seasonal winds are then averaged into 
different components to understand the role each component plays in the generation of 


eddies, jets and filaments. 


III. MODEL DESCRIPTION 


The numencal model used in this research was developed by Haney (1974,1985) 


for a closed basin, and modified by Batteen (1989) and Batteen e7 a/. (1989) for use in 


limited area eastern boundary current (EBC) regions with open borders on the northern, 


western and southern borders. It is explained in detail in Batteen e/ a/. (1989) and will be 


summarized below. 


A. MODEL EQUATIONS 


The model is a multi-level, primitive equation (PE) model, and uses the 


hydrostatic, ngid lid, B-plane and Boussinesq approximations. The governing equations 


are the following: 





a. Momentum Equations: 
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b. Continuity Equation: 
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d. Equation of State: 
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e: Thermodynamic Equation: 


dT , oui 
= Au VT + Ky 550, +87). (3.6) 


In the above equations, t is time, (x,y,z) 1S a mght-handed Cartesian coordinate system 
with x pointing eastward toward shore, y pointing northward alongshore, and z upward. 
The velocity components are (u,v,w), respectively. T is temperature, p is density and p’ is 
the departure from the vertically averaged pressure. & is a dummy vanable in equations 
(3.3) and (3.4). The density function is assumed to be a function of temperature only in 
equation (3.5). This is consistent with the region being modeled (Lynn ef al, 1982). 
Salinity may be a good tracer for water masses in the CCS (Huyer and Kosro, 1987); 
however it is not essential for a zero-order description since there are no major sources or 
sinks of salinity in this region of the CCS. In equation (3.6), O, = GS/p,Cc is the 


heating due to solar radiation, where 


S=S,(Re* +(1- R)e**). (3) 


S, is the downward flux of the solar radiation at the surface, R =.62 is the fraction of 
solar radiation absorbed in the upper few meters (z] = 1.5 m) and (1- R) =.38 is the 
fraction that penetrates to somewhat deeper levels (z2 = 20 m) as given by Paulson and 
Simpson (1977). The vertical turbulent mixing of heat and momentum by a dynamic 
adjustment mechanism is represented by the terms 6,(u), 6,(v) and 6,(7). This 
adjustment, a generalization of the convective adjustment mechanism, is based on the 
assumption of a cntical Richardson number, and it serves to maintain dynamic stability in 


the water column (Adamec e/ al., 1981). 


The boundary conditions at the top (z = 0)of the model are: 
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and at the bottom (z = - H) they are 
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t™ and tr are the cross-shore and alongshore components of the surface stress in 


equations (3.8a) and (3.8b). In (3.8c), Qp is the net upward flux of longwave radiation, 


sensible and latent heat across the sea surface. y = 10° is a geostrophic inflow angle in 


equations (2.9a) and (2.9b) (Weatherly, G. L., 1972). Table 1 provides constants used in 


this study as well as other symbols in the model equations. 


B. DOMAIN SIZE AND RESOLUTION 

The domain of the model is a rectangular region extending from 35°N to 45°N 
and 6° longitude off the west coast of the United States. This region extends 
approximately 500 km offshore and it spans over 1000 km in the alongshore direction. 
Point Conception marks the southern boundary of the study while Cape Mendocino is 
located in the middle of the study area. The honzontal resolution of the model is 8 km in 
the cross-shore direction and 17.5 km in the alongshore direction. This resolution should 
allow realistic spatial resolution of mesoscale features in the CCS, which have typical 
wavelengths of the order of 100 km (Breaker and Mooers, 1986). Variations in the 
coastline and ocean depth are omitted in this process-onented study to focus on the role of 


wind forcing in the generation of eddies and filaments. 


C. FINITE DIFFERENCE SCHEME AND TIME STEPPING 

A space staggered B-scheme which has u and v defined at one set of grd-points 
and T, p, and p at another (Arakawa and Lamb, 1977; Batteen and Han, 1981) is used in 
the horizontal. There are 10 layers in the vertical, separated by constant z-levels at depths 
of 13, 46, 98, 182, 316, 529, 870, 1416, 2283 and 3656 m. The time stepping used 
consists of a Matsuno time step followed by 10 leapfrog time steps and is continuously 


repeated throughout the model run. 


D. HEAT AND MOMENTUM DIFFUSION 
Biharmonic lateral momentum and heat diffusion with the coefficients listed in 


Table 1 are used in the model to allow the generation of mesoscale eddies. Laplacian 


lateral heat diffusion is not used since it can diminish the baroclinic signal associated with 


mesoscale processes (Holland and Batteen, 1986). 


E. SURFACE THERMAL FORCING 

The surface thermal forcing was highly simplified in this model in order to focus on 
wind's role in generating thermal variability. The solar radiation at the sea surface, So, 
was specified to be the mean summer CCS value from Nelson and Husby (1983). The 
mean summer CCS value was used even though the study used climatological winds for all 
Seasons because during the summer upwelling season there will be a net heat gain due to 
the cold upwelled water. Qp, the net longwave radiation was computed during model 
experiments from standard bulk formulas using the summer CCS mean value of 
alongshore wind, cloud cover, air temperature, relative humidity, and model-predicted sea 
surface temperature (Haney e/ a/, 1978). The initial sea-surface temperature was chosen 
so that the total heat flux across the sea surface, S, - Qp, equals zero. The only surface 
heat flux forcing was developed in Qp as a result of wind-forced fluctuations in the sea 
surface temperature. A more detailed discussion of this formulation can be found in 


Batteen ef al. (1989) and Haney (1985). 


F, HORIZONTAL BOUNDARY CONDITIONS 

The California coastline, the eastern boundary, is modeled as a straight vertical 
wall. A no-slip condition is imposed on the tangential velocity. 

The other three boundanes are open using a modified version of the radiation 
boundary conditions of Camerlengo and O'Brien (1980). Forcing to the boundaries 
caused unrealistic surface currents with no undercurrent in cases where the wind does not 


have an alongshore component. To remedy this problem, we impose a band where there 


is no wind forcing (McCreary, 1981; Batteen ef a/., 1989; McCreary ef a/., 1987). The 
bands of no wind forcing are located approximately 100 km from the northern and 
southern boundaries. This allows for the propagation of coastal Kelvin waves. A surface- 
trapped coastal jet and a relatively realistic undercurrent are established by the alongshore 


pressure gradient field generated by the Kelvin waves (Batteen ef al/., 1989). 


G. INITIAL CONDITIONS 
An exponential temperature profile with a vertical length scale of h = 450 m was 


the initial mean stratification. The form was: 


T =7,+VIe™, (3.10) 


where Tp = 2°C is the temperature at the bottom layer and AT = 13°C 1s the increase in 
temperature between the bottom and surface layer. This profile 1s considered to be 


climatologically representative of the CCS (Blumberg and Mellor, 1987). 


H. WIND DATA DESCRIPTION 

The winds used to force the model are from the European Centre for Medium 
Range Weather Forecasts (ECMWF) (Trenberth e/ a/., 1990) and were provided by the 
National Center for Atmosphenc Research (NCAR). World-wide mean monthly wind 
velocities were provided at 2.5 degree spacing for the years 1980 to 1989. Two 
climatologies denved from 80 months and 120 months of data were also provided. Since 
a comparison of the two climatologies showed no apparent contrast 1n our area of interest, 
the 120-month climotology was chosen to be the representative winds for this study. The 
data was linearly interpolated in both time and space to provide daily forcing at the 


required gnd points. The interpolated values are separated by approximately 8 km in the 
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cross-shore direction and 17.5 km in the along-shore direction, consistent with the model 
honzontal grid resolution. All the seasonal expenments used this 120-month climatolgical 


wind as the base for which the averaging (explained below) was done. 


I. EXPERIMENTAL DESIGN 
The first experiment is forced with a steady -830 cm s-! southward wind. The 
subsequent three time-varying cases broke the wind fields into different components to 
identify the importance each of these components played in the development of eddies and 
filaments. For the second and third experiments, the cross-shore component of the wind 
field 1s set to zero. The alongshore component in the second and third expenment ts then 
averaged and applied to the model gnd 1n the following manner: 
e varying in time but averaged spatially over the entire domain (second expenment); 
and 
° varying in time and in the alongshore direction but averaged in the cross-shore 
direction (third expenment). 
In the final experiment, the full climatological wind field is used. All four experiments are 


run for 360 days. 


J. ENERGY ANALYSIS TECHNIQUE 

The energy technique used in Batteen and Rutherford (1990) and Batteen ef al. 
(1992) will be used to analyze the generation and stability of the California Current 
System, eddies and filaments. The following is a summary of the Batteen et a/. (1992) 


descnption of this technique. 
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To gain an understanding of the energy transfers in the unstable flow, an energy 
analysis based on that of Han (1975) and Semtner and Mintz (1977) is made. The energy 
calculations with Semtner and Mintz (1977) notation 1s used in the energy calculations: 

QO time average 
() time deviation 
(") horizontal space average 
()° horizontal space deviation 
The kinetic energy 1s calculated by: 
(3.11) 


The time mean and time eddy kinetic energy are calculated, after reaching a quasi-steady 





State, by: 
—2 ,=2 
je ae (3.12) 
2 
| aes (3.13) 
2 
Available potential energy (P) ts: 
8 
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This determines when a quasi-steady state 1s reached in order to collect statistics. The 


temporal mean and eddy available potential energy are then calculated by: 


P = oe Cie uy (3.15) 


P! = oF (Ty TFS (3.16) 


Following Semtner and Mintz (1977), transfers between energy types are defined as: 


{K + P}=-og|7w] G19) 





|P'> K'}= ag|T'w"| (3.18) 
[R>K'}=0(ver¥ +0] (3.19) 
{P— P'}=ag| T’Vew'T’ {Z| , (3.20) 


Vertical velocity and advection terms calculated by the model are not stored. Vertical 
velocity and advection terms are recalculated in the same manner as in the model using 
stored values of u, v and T as the input data. The recalculated energy transfers are 
consistent with the initial calculations of vertical velocity and advection terms in the 
model. 

The energy transfer analysis is used to argue for the instability mechanism which 
leads to the initial eddy generation by studying the quasi-steady energetic state prior to 


eddy generation. 


Pa 


Po 


Table 3.1. VALUES OF CONSTANTS USED IN THE MODEL 


VALUE 


0.958 cal gm™! @Ky! 
1.225x 10° 
278.2°K 


1.23 x 10° pm cm 


1.0276 pm cm? 
2.01 x 104 (Ky! 
10 

7.6 x 10° cm 

1.75 x 10° cm 
4.5x 10° cm 
800s 

0.93x 104s" 
980 cms’! 

2x 10!7 om4 s7] 
2x 101? cm4 s7) 
0:5 cm s! 


0.5 cm? s] 


DEFINITION 


specific heat of sea water 

bottom drag coefficient 

constant reference temperature 
density of air 

density of sea water at Ty 

thermal expansion coefficient 
number of levels in vertical 
cross-shore grid spacing 

alongshore gnd spacing 

total ocean depth 

time step 

mean Coriolis parameter 
acceleration of gravity 

biharmonic momentum diffusion coefficient 
biharmonic heat diffusion coefficient 
vertical eddy viscosity 


vertical eddy conductivity 
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IV. RESULTS OF WIND FORCING EXPERIMENTS 


Experiment | studies the effects of steady, equatorward (upwelling-favorable) 
winds in the CCS. Experiments 2 through 4 investigate the effects of seasonally varying 
winds in the same region, using climatological monthly values that have been interpolated 
to daily values. Experiment 2 uses temporally varying, spatially averaged winds, while 
Experiment 3 uses winds with temporal vanation and alongshore vanability. In 


Experiment 4 the full climatological, seasonally varying wind field is used. 


A. EXPERIMENT 1 

In Experiment 1] the model was forced with a steady -830 cm s-! equatorward 
wind. This wind forcing of 1 dyne cm-2 is the summer (June-August) average for the 
modeled area (Batteen ef al, 1989). Inertial oscillations of near-surface ocean currents 
developed initially, as in Batteen ef al. (1991). After a few days the oscillations were 
damped, leaving quasi-steady offshore Ekman transport (not shown). By day 20 an 
equatorward coastal jet developed with speeds of ~15 cm s7! (Fig. 4.la,c). Colder 
surface temperatures also appeared near the coast (Fig. 4.1b). A poleward undercurrent 1s 
also discernible at 182 m depth (Fig. 4.1d). The surface current flowed at ~25 cm s"!\ and 
extended ~35 km offshore (Fig. 4.1c), while the undercurrent, centered along the coast, 
had weaker velocities of ~5 cm s~! (Fig. 4.1d). 

Plots of the zonal velocity (Figs. 4.2a, b) show that eddies begin to form around 
days 60 and 70 as evidenced by perturbations in the zonal velocity fields. (Since 
meridional velocity plots tend to be dominated by the predominately meridional currents, 
perturbations are more clearly seen in zonal rather than mendional velocity plots.) The 


isotherms also show meanderings in the same region where the eddies are generated 
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(Figs. 4.2c, d). As expected, eddies first formed in the poleward end of the region, as in 
the B-plane expenment of Batteen ef a/. (1989). The vertical cross-section of mendional 
velocity at day 90 (Fig. 4.3) shows that the surface equatorward current has its core at 
~80 m depth with a speed of ~25 cm s-!,, while the undercurrent develops below a depth 
of ~240 m within ~25 km of the coast. 

Both baroclinic and barotropic transfers (equations 3.18 and 3.19, respectively) 
can contnbute to the formation of eddies. Baroclinic instability draws its energy from the 
vertical shear in the mean flow while barotropic instability receives energy from the 
horizontal shear (Haltiner and Williams, 1980; Batteen ef a/, 1989). Figure 4.4 shows that 
although both baroclinic (Fig. 4.4a) and barotropic (Fig. 4.4b) instabilities play a key role 
in the generation of eddies, barotropic instability 1s the dominant generation mechanism. 
Plots of the surface velocity vectors supermmposed on the isotherms (Fig. 4.5a-d)and the 
menidional velocity contours (Fig. 4.6a-d) at 90 day intervals show how the meandenng 
current continues to grow in strength and to spin off isolated eddies. This strong 
meandering current can produce eddies that have wavelengths up to ~ 200 km and that 
can propagate at least ~ 200 km offshore (Figs. 4.5 and 4.6). 

Although these results show a meandering current and eddies, it 1s very unlikely 
that there will be steady strong equatorward (upwelling favorable) winds for 360 days, 
since there are many penods of poleward winds and also wind relaxations in the CCS 
region. Therefore, three more experiments with a time varying wind, were conducted in 
order to better understand what cntical unsteady wind forcing elements are needed for the 


formation and maintenance of currents and eddies. 


18 


B. EXPERIMENT 2 

The second expenment utilized the alongshore-component of the temporally 
averaged ECMWF winds. The seasonal cycle of this wind field is plotted in Figs. 4.7 and 
4.8. At the beginning of the year, the wind forcing was very weak, and the winds 
remained weak until day 135. From days 135 to 225, the wind was equatorward and 
continued to grow in strength. The wind then decreased in strength until it was nearly 
zero again at day 360. 

The model simulation results show an equatorward coastal jet and colder surface 
temperatures near the coast by day 180 (Fig. 4.9b). The surface current flowed at ~10 cm 
s-! and extended ~25 km offshore, while the undercurrent had weaker velocities of ~1.5 
cm s~! and was centered along the coast at ~120 m depth (Fig. 4.10). 

The isotherms show meanderings in the same region where the eddies are 
generated (Figs. 4.9 and 4.10). Again the eddies first formed in the poleward end of the 
region, aS expected. Fig. 4.12 shows that both baroclinic and barotropic instabilities 
contribute to the generation of eddies, with the barotropic instability being the dominant 
mechanism. The plots of the surface velocity vectors supernmposed on the isotherms at 
selected days between 130 and 270 show how the meandemnng current continues to grow 
in strength and to spin off isolated eddies (Figs. 4.9 and 4.10). As expected, this 
meandering current weakens and dissipates as the seasonal cycle of the temporal wind 
weakens. 

With this non-spatially varying seasonal cycle, currents and eddies are generated 
much later than with the ideal upwelling favorable winds of the first experiment. The 


major portion of activity happened between day 135 and day 290 and only showed very 
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weak upwelling and little eddy activity (Figs. 4.9-4.10). The eddies also only propagated 
to ~ 100 km offshore. 


C. EXPERIMENT 3 

The third experiment uses the alongshore component of the time-dependent wind 
forcing with spatial vanability in latitude The wind forcing is divergent at the beginning of 
the year with poleward winds north of 40° N latitude and equatorward winds south of 40° 
N latitude. The wind stress shifts to being equatorward everywhere by day 135, 
intensifying from days 135 to 225, and then splits again by day 315 (Figs. 4.13 and 4.14). 

The ocean response 1s evident by day 70 when an equatorward coastal jet forms 
(not shown). Eddies form by days 240 and 260 (Figs. 4.15a,b). Meanders in the 
isotherms in the areas of these eddies are also evident (Figs. 4.15c,d). The cross-section 
of the meridional velocity at day 260 (Fig. 4.16) shows an equatorward surface current 
with speeds of ~35 cm s"! and a poleward undercurrent with its core located 40 km 
offshore and at a depth of 350 m. 

In this experiment, there is a competition between the B-plane effect, which will 
cause eddies to form in the poleward end of the domain, and the continuous and stronger 
equatorward wind forcing in the southern portion of the domain. This causes eddies to 
form throughout the region so that by day 300 the eddies encompass the area from the 
northern boundary to the southern boundary of the domain (Figs. 4.17a,c). Day 360 
shows the strong meandering current with large eddies imbedded in it (Figs. 4.17b,d). 

Figures 4.18a,b show that while both baroclinic and barotropic instabilities are 
eddy generating mechanisms in the poleward portion of the domain, it 1s primarily the 
barotropic instability that contnbutes to the generation of eddies in the southern portion of 


the domain where the wind is predominately equatorward. The wind forcing in this 
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experiment produced larger eddies (that propagate at least ~ 150 km from shore) and eddy 
velocities than were seen in the previous expemment, which is more consistent with 


available observations. 


D. EXPERIMENT 4 

The final experiment uses full climatological wind forcing (Figs. 4.19 and 4.20). 
Again the divergence in the winds at 40° N is evident in the beginning of the year. By day 
135, the winds over the entire domain have an equatorward component and by day 225 
the Pinds are predominately equatorward everywhere and strongest in magnitude. After 
day 225 the winds become less intense and by day 360, north of 40° N the wind becomes 
poleward again. 

As in the previous cases, inertial oscillations of near-surface ocean currents 
developed initially. After a few days these observations were damped, leaving quasi- 
steady offshore Ekman transport (not shown). By day 20 an equatorward coastal jet 
developed(not shown). At day 110 the surface current flowed at ~40 cm s"!, and 
extended ~35 km offshore, while a weak undercurrent was evident (Fig. 4.21). The 
meridional velocity cross-section at day 180 (Fig. 4.22) shows a strong equatorward 
surface jet centered ~50 km offshore with speeds of ~60 cm s7!. The poleward 
undercurrent is located at ~360 m depth and ~40 km offshore. 

Plots of the zonal velocity (Fig. 4.23d) show that eddies begin to form around day 
260 as evidenced by perturbations in the zonal velocity fields. As expected, the isotherms 
also show meanderings in the same region where the eddies are generated (Fig. 4.23c). 
The eddies again formed initially in the poleward end of the region. 

Figure 4.24 shows that, although both baroclinic (Fig. 4.24a) and barotropic (Fig. 


4.24b) instabilities play a key role in the generation of eddies, barotropic instability is the 
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dominant generation mechanism. Days 300 to 360 (Figs. 4.25-4.27) show the current 
becoming stronger with the eddies continuing to grow and to separate from the coast. 
The speeds of the currents and eddies reach ~ 65-85 cm s“! with wavelengths of ~ 300 km 
and can propagate to ~ 350 km offshore. This case gives results that are more consistent 
with available observations than was obtained in any of the previous cases. However, the 


speed and depth of the undercurrent are still not quite in agreement with observations. 
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Figure 4.1a) Experiment 1 Day 20: Surface velocity vectors superimposed on isotherms. The contour 
interval for the isotherms is 0.59 C. The temperature decreases toward the coast. 
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Temperature conlour at model day 20.0 


Figure 4.1 b) Experiment 1: Day 20 Temperatures. The contour interval for the isotherms is 1.0° C. 
The temperature decreases toward the coast 
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Figure 4.1 c) Experiment 1: Day 20 Meridional velocities. The contour interval is 10 cm s"?. The dashed 
lines indicate southward velocities. The equatorward velocity increases to ~ 25 cm s”! near the coast. 
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COKTOUR FROM ~11 TO 10 BY 3 
North-South vel contour at model day 20.0 


Figure 4.1 d) Experiment 1: Day 20 Meridional velocites at 182m. = The contour interval is 3 cm s“*. 
The dashed lines indicate southward velociues. The poleward undercurrent 1s seen at the coast. 
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East-West vel contour at model day 60.0 


Figure 4.2 a) Experiment 1: Surface Zonal velocities at day 60. The contour interval is 10 cm s"*. The 
dashed lines indicate westward velocities. Closed contours show the locations of eddies. 
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East-West vel contour at model dey 70.0 


Figure 4.2 b) Expermment 1: Surface Zonal velocities at day 70. The contour interval is 10cm s**. The 
dashed lines indicate westward velociues. Closed contours show the locauons of eddies. 
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Figure 4.2 c) Experiment 1: Surface temperatures at day 60. The contour interval is 19 C. The 
temperature decreases toward the coast. 
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Figure 4.2 d) Experiment 1: Surface temperatures at day 70. The contour interval is 1° C. The 
temperature decreases toward the coast. 
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Figure 4.3 Experiment 1: Cross-section at 40° N latitude of meridional velocities at day 90. The contour 
interval is 1.1 cm s"!. The dashed lines indicate southward velocities and show the surface equatorward 
current. The solid lines indicate northward flow and show the poleward undercurrent. 
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Figure 4.4 Experiment 1 Transfers of energy (as defined in equations 3.18 and 3.19) from: a) P toK 


(eddy available to eddy kinetic energy), and b) KtoK (mean to eddy kinetic energy). Transfers of 
energy are averaged over the days 50-70 and summed over the upper five layers. The contour intervals 
are 0.00005 ergs cm-3 s"! in a) and 50 ergs cm73 5°! in b). 
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Figure 4.5a) Experiment 1: Surface velocity vectors superimposed on temperatures at day 90. The contour 
interval is 0.5° C. The temperature decreases toward the coast. 
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Temp and velocity vector al model day 180.0 
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Figure 4.5 b) Experiment 1: Surface velocity vectors superimposed on temperatures at day 180. The 
contour interval is 0.59 C. The temperature decreases toward the coast. 


34 


Temp and velocity vector at model day 270.0 
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Figure 4.5 c) Experiment 1: Surface velocity vectors superimposed on temperatures at day 270. The 
contour interval is 0.59 C. The temperature decreases toward the coast. 
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Temp and velocily vector at model day 360.0 
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Figure 4.5 d)Expceriment 1: Surface velocity vectors superimposed on temperatures at day 360. The 
contour interval is 0.59 C. The temperature decreases toward the coast. 
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CONTOUR FROM -36 TO 38 BY 10 
North-South vel contour at model day 90.0 


Figure 4.6 a ) Experiment 1: Surface meridional velocities at day 90. The contour interval is 10 cm s*?. 
The dashed lines indicate southward velocities. Closed contours show the locations of eddies. 
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CONTOUR FROM -36 TO 56 BY 10 
North-South vel contour at model day 160.0 


Figure 4.6 b ) Experiment 1: Surface meridional velocities at day 180. The contour interval is 10 cm s“?. 
The dashed lines indicate southward velocities. Closed contours show the locations of eddies. 
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North-South vel contour at model day 270.0 


Figure 4.6 c) Experiment 1: Surface meridional velocities at day 270. The contour interval is 10 cm s"?. 
The dashed lines indicate southward velocities. Closed contours show the locations of eddies. 
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North-South vel contour at model day 360.0 


Figure 4.6 d) Experiment 1: Surface meridional velocities at day 360. The contour interval is 10 cm s7?. 
The dashed lines indicate southward velocities. Closed contours show the locations of eddies. 
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Figure 4.7a) Experiment 2: Wind forcing in m s~* at day 45. 
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Figure 4.7 b) Experiment 2: Wind forcing in m s~* at day 90. 
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Figure 4.7 c) Experiment 2: Wind forcing in m s~* at day 135. 
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Figure 4.7 d) Experiment 2: Wind forcing in m $s” 
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Figure 4.8 b) Experiment 2: Wind forcing in m s™* at day 270. 
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Figure 4.8 c) Experiment 2: Wind forcing in m s~* at day 315. 
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Figure 4.8 d) Experiment 2: Wind forcing in m s~* at day 360. 
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Temp and velocity vector at model day 130.0 
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Figure 4.9 a) Experiment 2: Velocity vectors superimposed on surface temperatures at day 130. The 
contour interval is 0.59 C. The temperature decreases toward the coast. 
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Temp and velocity vector at model day 180.0 
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Figure 4.9 b) Experiment 2: Velocity vectors superimposed on surface temperatures at day 180. The 
contour interval is 0.5° C. The temperature decreases toward the coast. 
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Temp and velocity vector at model day 2000 
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Figure 4.9 c) Experiment 2: Velocity vectors superimposed on surface temperatures at day 200. The 
contour interval is 0.5°C. The temperature decreases toward the coast. 
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Temp and velocity vector at model day 210.0 
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Figure 4.9 d) Experiment 2: Velocity vectors superimposed on surface temperatures at day 210. The 
contour interval is 0.5°C. The temperature decreases toward the coast. 


52 


1OSze 


87S 


TA TAy 


Latitude (km) 
ul 
a 
a) 


4@2. 


245. 


a7. 


Temp and velocity vector at model day 250.0 
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Figure 4.10 a) Experiment 2: Velocity vectors superimposed on surface temperatures at day 250. The 
contour interval is 0.59 C. The temperature decreases toward the coast. 
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Temp and velocity veclor at model day 270.0 
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Figure 4.10 b) Experiment 2: Velocity vectors superimposed on surface temperatures at day 270. The 
contour interval is 0.5°C. The temperature decreases toward the coast. 
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Temp and velocity vector al model day 280.0 
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Figure 4.10 c) Experiment 2: Velocity vectors superimposed on surface temperatures at day 280. The 
contour interval is 0.5° C. The temperature decreases toward the coast. 
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Temp and velocity vector at model day 290.0 


1832. 


8752 


PAT 


é 

= S6Ge 

@ 

° 

= 

s 80 cm/s 
MAXTHUM VECTOR 


4Be. 


245. 


87. 





-200.80 -166.7 -133.3 -108.@ -66.7 -33.3 9.2 
Distance off shore (km) 


Figure 4.10 d) Experiment 2: Velocity vectors superimposed on surface temperatures at day 290. The 
contour interval is 0.59 C. The temperature decreases toward the coast. 
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Figure 4.11 Experiment 2: Cross-section at 40° N latitude of meridional velocities at day 210. The 
contour interval is 1.1 cm sl, The dashed lines indicate southward velocities and show the surface 
equatorward current. The solid lines indicate northward flow and show the poleward undercurrent. 
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Figure 4.12 Experiment 2 Transfers of energy (as defined in equations 3.18 and 3.19) from: a) P wok 


(eddy available to eddy kinetic energy), and b) KtoK (mean to eddy kinetic energy). Transfers of 
energy are averaged over the days 300-320 and summed over the upper five layers. The contour intervals 
are 0.0000025 ergs cm*> s"! in a) and 10 ergs cm”? s“! in b). 
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Figure 4.13a) Experiment 3: Wind forcing in m s~* at day 45. 
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Figure 4.13 b) Experiment 3: Wind forcing in ms~* at day 90. 
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Figure 4.13 c) Experiment 3: Wind forcing in m s~* at day 135. 
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Figure 4.13 d) Experiment 3: Wind forcing inm s” 


at day 180. 


Figure 4.14 a) Experiment 3: Wind forcing in m s” 
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Figure 4.14 b) Experiment 3: Wind forcing in m s~* at day 270. 
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Figure 4.14 c) Experiment 3: Wind forcing in m s! at day 315. 
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Figure 4.14 d) Experiment 3: Wind forcing in m s"* at day 360. 
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CONTOUR FROM -26 TO 88 BY 10 
East-West vel contour at model day 240.0 


Figure 4.15 a) Experiment 3: Surface Zonal velocities at day 240. The contour interval is 10 cm s~’*. The 
dashed lines indicate westward velocities. Closed contours show the locations of eddies. 
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Figure 4.15 b) Experiment 3: Surface Zonal velocities at day 260. The contour interval is 10 cm s“°. 
The dashed lines indicate westward velocities. Closed contours show the locauons of eddies. 
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CONTOUR FROM 10 TO 16 BY 1 


Temperature contour at model day 240.0 


Figure 4.15 c) Experiment 3: Surface temperatures at day 240. The contour interval is 19 C. The 
temperature decreases toward the coast. 
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Figure 4.15 d) Experiment 3: Surface temperatures at day 260. The contour interval is 19 C. The 
temperature decreases toward the coast. 
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North-South velocity at model day 260.0 
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Figure 4.16 Experiment 3: Cross-section at 40° N latitude of meridional velocities at day 260. The 
contour interval is 1.1 cm s*!. The dashed lines indicate southward velocities .and show the surface 
equatorward current. The solid lines indicate northward flow and show the poleward undercurrent. 
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Temp end velocity veclor al model day 3000 
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Figure 4.17 a) Experiment 3: Surface velocity vectors superimposed on temperatures at day 300. The 
contour interval is 0.59 C. The temperature decreases toward the coast. 
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Figure 4.17 b) Experiment 3: Surface velocity vectors superimposed on temperatures at day 360. The 
contour interval is 0.59 C. The temperature decreases toward the coast. 
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Figure 4.17 c) Experiment 3: Surface temperatures at day 300. The contour interval is 0.59 C. The 
temperature decreases toward the coast. 
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Figure 4.17 d) Experiment 3: Surface temperatures at day 360. The contour interval is 0.5° C. The 
temperature decreases toward the coast. 
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Figure 4.18 Experiment 3 Transfers of energy (as defined in equations 3.18 and 3.19) from: a) P to K 


(eddy available to eddy kinetic energy), and b) KtoK (mean to eddy kinetic energy). Transfers of 
energy are averaged over the days 280-310 and summed over the upper five layers. The contour intervals 
are 0.00025 ergs cm} s-! in a) and 1000 ergs cm73 5"! in b). 
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Figure 4.19a) Experiment 4: Wind forcing in m s~* at day 45. 
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Figure 4.19 b) Experiment 4: Wind forcing in m s~* at day 90. 
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Figure 4.19 c) Experiment 4: Wind forcing inm s“* at day 135. 
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Figure 4.19 d) Experiment 4: Wind forcing inm $s" 


at day 180. 


— Le Le LO en ees eee Se —— —  —— —— —— 
i Oc Wc cD Pee De ee Dee eee ae ee oe oe oe 
Ww ew UD Ue Oe Oo ee? ee? ae ee oo = — > 
TF Fee Oe Oe ee > ee >  —— ==> -—-—— > 
a Eg ee ee Ep eee ees Se ee 


Figure 4.20 a) Experiment 4: Wind forcing in m s” 


at day 225. 
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Figure 4.20 b) Experiment 4: Wind forcing in m s~* at day 270. 
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Figure 4.20 c) Experiment 4: Wind forcing in m s~* at day 315. 
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Figure 4.20 d) Experiment 4: Wind forcing in m s"* at day 360. 
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Figure 4.21 a) Experiment 4: Surface velocity vectors superimposed on temperatures at day 110. The 
contour interval is 0.5° C. The temperature decreases toward the coast. 
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Figure 4.21 b) Experiment 4: Cross-section at 40° N latitude of meridional velocities at day 110. The 
contour interval is 1.1 cm s"!. The dashed lines indicate southward velocities and show the surface 
equatorward current. The solid lines indicate northward flow and show the poleward undercurrent. 
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Figure 4.22 Experiment 4: Cross-section at 409 N latitude of meridional velocities at day 180. The 
contour interval is 1.1 cm s7}. ‘The dashed lines indicate southward velocities and show the surface 
equatorward current. The solid lines indicate northward flow and show the poleward undercurrent. 
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North-South vel contour at model day 260.0 


Figure 4.23 b) Experiment 4: Day 260 Meridional velocities. The contour interval is 10 cm s™*. The 


dashed lines indicate southward velocities. The equatorward velocity increases to ~35 cm s"! near the 
coast. 
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Temperature contour al model day 260.0 


Figure 4.23 c) Experiment 4: Surface temperatures at day 260. The contour interval is 19 C. The 
temperature decreases toward the coast. 


90 


Distance off shore (km) 
-512.8 -426.7 -341.3 -256.8 -178.7 -865.3 @.@ 


1832.5 


675.8 


717.5 


562.2 


Latitude (kam) 


402.5 


245.8 


67.5 





CONTOUR FROM -25 TO 8S BY 10 
East-West vel contour at model day 260.0 


Figure 4.23 d) Experiment 4: Surface Zonal velocities at day 260. The contour interval is 10 cm s"+. The 
dashed lines indicate westward velocities. Closed contours show the locations of eddies. 
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Figure 4.24 Experiment 4 Transfers of energy (as defined in equations 3.18 and 3.19) from: a) P toK 


(eddy available to eddy kinetic energy), and b) Kwok (mean to eddy kinetic energy). Transfers of 
energy are averaged over the days 230-250 and summed over the upper five layers. The contour intervals 
are 0.0005 ergs em-3 57! in a) and 2000 ergs cm} 5°! in b). 
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Figure 4.25 a) Experiment 4: Surface velocity vectors superimposed on temperatures at day 300. The 
contour interval is 0.5° C. The temperature decreases toward the coast. 
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Figure 4.25 b) Experiment 4: Day 300 Meridional velocities. The contour interval is 20 cm s™*. The 
dashed lines indicate southward velociues. 
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Figure 4.25 c) Experiment 
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Figure 4.25 d) Experiment 4: Surface Zonal velocities at day 300. The contour interval is 10cm s"*. The 
dashed lines indicate westward velocities. Closed contours show the locauons of eddies. 
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Figure 4.26 a) Experiment 4: Surface velocity vectors superimposed on temperatures at day 340. The 
contour interval is 0.59 C. The temperature decreases toward the coast. 
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North-South vel contour at model day 340.0 


Figure 4.26 b) Experiment 4: Day 340 Meridional velociues. The contour interval is 20 cm s**. The 


dashed lines indicate southward velociues. 
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Figure 4.26 c) Experiment 4: Surface temperatures at day 340. The contour interval is 1° C. The 
temperature decreases toward the coast. 
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CONTOUR FROW -28 TO 26 BY 10 
East-West vel contour at model day 340.0 


Figure 4.26 d) Experiment 4: Surface Zonal velocities at day 340. The contour interval is 10 cm s~*. The 
dashed lines indicate westward velocities. Closed contours show the locations of eddies. 
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Figure 4.27 a) Experiment 4: Surface velocity vectors superimposed on temperatures at day 360. The 
contour interval is 0.59 C. The temperature decreases toward the coast. 
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CONTOUR FROM -156 TO 85 BY £0 
North-South vel contour at model day 360.0 


Figure 4.27 b) Experiment 4: Day 360 Meridional velocities. The contour interval is 20 cm s“*. The 
dashed lines indicate southward velociues. 
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Figure 4.27 c) Experiment 4: Surface temperatures at day 360. The contour interval is 19 C. The 
temperature decreases toward the coast 


103 


Distance off shore (km) 
-512.@ -426.7 -341.3 -256.8 -178.7 -85.3 0.2 


1832.5 





675.8 


719% 


Latitude (km) 
ul 
oO 
& 
a 


402.5 


245.2 


67.5 


CONTOUR FROM -86 TO 85 BY 10 
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Figure 4.27 d) Experiment 4: Surface Zonal velocities at day 360. The contour interval is 10 cm s“?, 
The dashed lines indicate westward velocities. Closed contours show the locauons of eddies. 
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V. SUMMARY AND RECOMMENDATIONS 


A. SUMMARY 

In all four cases, currents and eddies were generated which shows the significant 
role wind forcing plays in the generation of these features. Table 5.1 shows a summary 
comparing the model results with available observations. The results of Experiment 1 are 
fair except for the weak surface current and undercurrent. Though Experiment 1 shows 
mostly reasonable results compared with observations, it is not likely that there would be 
360 days of strong equatorward (upwelling favorable) winds. Experiment 2 showed the 
least reasonable results compared with observations. The coastal jet and undercurrent 
were much weaker and smaller than observations in Expenment 2. In all four 
experiments, the undercurrent was wider than seen in observations. A possible reason the 
undercurrent is wider than seen in observations may be due to not having topography 
representing the continental shelf in the model. Experiment 3 showed reasonable results 
with the exception of a very wide undercurrent and small eddy velocities. The results that 
were closest to observations were from Experiment 4. With the full climatological winds, 
the model was able to produce a coastal jet and eddy velocities of 85 cm s-!. The 
conclusion obtained from the results of the four experiments is that both temporal and 
spatial variability in wind forcing in eastern boundary current regions are critical elements 


in the formation and maintenance of currents and eddies. 


B. RECOMMENDATIONS 
Future studies should look at the role of inter-annual vanability in the wind 
forcing. By allowing the model to run for ten years (the 1980-1989 data set is available), 


it may be possible to detect what roles anomalous years such as E/ Nino have on the 
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currents and eddies in the CCS. Longer model simulations would also help to determine 
the life span of the larger eddies. 

Experiments to include the topography and an irregular coastline to study the role 
of these features in the maintenance of currents and eddies also may make it possible to 
obtain a more reasonable undercurrent. Nevertheless, since both currents and eddies have 
been generated and maintained in these wind forcing experiments, we must conclude that 
wind forcing is a critical element in the formation and maintenance of these features in the 
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Table 5.1 INSTANTANEOUS COMPARISON OF EXPERIMENTS (EXP.) WITH OBSERVATIONS 
(OBS.) OF THE CCS. 


_ Maximum coastal jet velocity (cm s"!) 


. Offshore location of coastal jet (km) 

. Offshore extent of coastal jet (km) 

. Depth of inshore coastal jet (m) 

. Maximum undercurrent velocity (cm sl) 


. Offshore location of undercurrent axis (km) 


Maximum width of undercurrent (km) 
Depth of undercurrent axis (m) 
Maximum zonal eddy diameter (km) 


Maximum zonal eddy velocity (cm $7!) 


(3) Kosro (1987) 
(4) Davis (1985) 


(1) Kosro and Huyer (1986) 
(2) Huyer and Kosro (1987) 


30-100 (1,2,3,4) 
25-35 (2.3) 
>40 (1,2,3) 
90-150 (2,3) 

5-15 (2,3) 
10-40 (2,3) 
10-20 (2,3) 
200-300 (2) 

10->100 (2,5,6,7) 


50-100 (1,2,3,4.6,7) 


(5) Mooers and Robinson (1984) 
(6) Brink and Cowles (1991) 


(7) Brink et al. (1991) 
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Exp. 
1 
25 
25-35 
100 
80 
3 
25 
20 
280 
145 


25 


Exp. 
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10-15 


20 


60 
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20 


25 


100 


40 


Exp. 


35 
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40 
80 
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40 
80 
350 
120 


35 
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